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1 Introduction 

^ We present a uniform reconstruction method for a class of Funk-Radon type integral 

transforms extending results of [16] for arbitrary dimensions. The reconstruction does 
^ not include summation of an infinite series and looks as the standard (John type) inversion 

\Q for the Radon transform. We specify this method for classical and new regular acquisition 

, J geometries. The condition of regularity (see the next Section) is necessary for an inversion 

^ operator to be bounded in a Sobolev space scale, but it is not sufficient. We shall see 

I that the existence of a elementary exact reconstruction formula depends on vanishing 

of some integrals of rational forms on an algebraic manifold. In the last Section, we 
discuss reconstruction in Euclidean space from data of integrals over a family of spheres. 
This subject is in focus of recent research, see papers [1] , [10] , [H] , [19] , [I3j and references 



X 

^ therein. 



2 Geometry and integrals 

Let X and S be smooth n dimensional manifolds (ra > 1), let Z be a smooth closed 
hypersurface in A x S and p : Z X, 7r:Z— J-Ebe natural projections. We suppose 
that there exists a real smooth function $ in A x S (called generating function) such 
that Z = {{x,(t) ; $ (x, a) = 0} and d^^ 7^ on Z. Suppose that 

(i) The map it has rank n and the mapping P : N* (Z) T* (A) is a local diffeomor- 
phism. Here, N* (Z) denotes the conormal bundle of Z and P (x, a; z/j,, z/g.) = {x,i>x) e 
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T* (X) . It follows that a set Z io) = (a) = {x; $ (x, a) = 0} is a smooth hypersur- 
face in X, and for any point x G X and for any tangent hyperplane h G {X) there is 
a locally unique hypersurface Z (A, u) through x tangent to h. 

Proposition 2.1 For an arbitrary generating function $ property (i) is equivalent to the 
condition: det (da;^jda-,r^) 7^ where \& (x, t; a, r) = tr^ (x, cr) , t, r G M, tr > /or any 
local coordinate system X and any local coordinate system ai, m S. 

For a proof see [I5], Proposition 1.1. 

Definition. We call a generating function $ regular if it satisfies conditions (i) and 
(ii) there are no conjugated points, that is the equations $ (x, cr) = $ {y, cr) and 

do-$ (x, 0") = do-$ (y, 0") are fulfilled for no x 7^ y G X, a G S. 

Suppose that X is oriented and g is a Riemannian metric in X; let dV^ be the oriented 

Riemannian volume form. Consider the integral 



for an arbitrary continuous function / compactly supported in X. The quotient fdV/ d<I> 
denotes an arbitrary n — 1 form q such that d$ A q = fdV. It is defined up to a term 
hd6 where /i is a continuous function. An orientation of a hypersurface Z (a) is defined 
by means of the form d6 and the integral of q over Z (a) is uniquely defined. We call the 
operator Mq, Minkowski-Funk transform generated by $. This transform can be written 
in terms of Riemannian integral as follows: 



Here, dgS is the Riemannian n—1 surface element and |Vga| = y g (da) is the Riemannian 
gradient of a function a. The function |Vg$| vanishes nowhere since of (i). Suppose that 
the gradient factorizes through X and S that is |Vg$(x,cr)| = m{x)fi{a) for some 
positive continuous functions m in X and /i in S. Then data of the Minkowski-Funk 
transform is equivalent to data of Riemann hypersurface integrals 



since Rgf {a) = fi{(T) M,^ (mf) {a). The reconstruction problem of a function / from 
Riemann integrals Rgf is then reduced to inversion of the operator M$. 

We say that a generating function $ is resolved if S = M x S""^, and $ has the form 
$ (x; X,u) = 6 (x, u) — X, A G M, w G S"^"'^ for a smooth function 6* on X x S"^"'^ where 
S"'^^ denotes the unit sphere in an n dimensional space. Note that the map p : Z ^ X 
is always proper for a resolved generating function. This property guarantees that the 
functions M$/ and Rgf have compact support in E. The operator M$ fulfils the range 
conditions similar to that of the Radon transform. 





(1) 
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Proposition 2.2 Let ^ = 6 — \ be a resolved regular generating function and 9 {x,uj) be 
a polynomial function of cu of order m. Then for an arbitrary integrable function f in X 
with compact support and for an arbitrary polynomial p (A) of order k the integral 

p{\)M^f {\uj)d\ 



is a polynomial of u of order < mk. 
Proof. We have 



p{\)AUf(\,ui)d\= I p{X) I ^dA= / p(fl(i,w))/(x)dV 



=A JX 

where p {9 (x, uj)) is a polynomial of ui of order < mk. ► 

3 Main theorem 

For a real smooth function / in a manifold X and a natural n we consider singular 
integrals 

In± (p) = / . ./.^vn = lim ^ 



for a smooth density p with compact support. If d/ 7^ on the zero set of /, then these 
limits exist and the functionals /„± are generalized functions in X. The functional 



is called a principal value integral. For a resolved regular generating function ^ — 9 — \ 
we define the function on X x X\{diag} 

. . da; 



©n {x, y)^ 



g„-i {9 {x, (jj)-9 {y, uj) - iO)'' 

where do; is the Euclidean volume form on S"^~^. The singular integral converges since 
by (ii) the d^ {9 {x, uj) - 9 {y, u)) ^ as 9 {x,uj) - 9 {y, u) = 0. 

Theorem 3.1 Let ^ = 9 — X be a regular resolved generating function in X xJ] and f G 
L2 {X) be an arbitrary function with compact support. If n is even and Re 0„ {x, y) = 
for any x ^ y E X, a reconstruction from data of Mq>f is given by the formula: 

_ {n-l)\ 1 f M^f{\uj)dMu 

, , I (P) A 
{2m) Dn [x) Jsn-i Jk 9A" ^ 9{x,u)- X 
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If n is odd and Im0„ {x,y) = for x ^ y, the function can be reconstructed by 

f (^) = , . ^ , , I 5^""'^ {0 ix, uj) - A) M^f (A, uj) dAda; (3) 
2 27rz D„ ix) Jt. 



2{2tii 

1 f (9"-i 



2 (27rz)""^ Dn (x) Js"-! 
D„ fx) = 



-M$/(A,a;) 



do; 

A=e(a:,t^) 



dw 



n-i |Vg6' (x, a;)!"' 



The integrals and converge in mean on any compact set in X . 

Remark 1. A more invariant form of ([2]) or ([3]) is the reconstruction of the form 
fdV: 

{2in} Dn J 

The quotient dV/Dn depends only on the conformal class of the Riemannian metric g. 
Remark 2. The condition (ii) is not required when n = 2. 

Remark 3. G. Beylkin studied "the generalized Radon transform" [1], which coin- 
cides with the operator M$ in the Euclidean space. He constructed a Fourier integral 
operator parametrix for this operator and reduced inversion of this operator to solution 
of a Fredholm equation. Two dimensional case of Theorem 3.1 was stated in [16j. 

Remark 4. When the condition on 0„ fails, the right-hand side Rn of ([2]), respec- 
tively Q, gives anyway a reconstruction up to a compact operator. More precisely, let 
K be an arbitrary compact set in X and Ik be the indicator function of K. Then the 
equation holds IxRf = / + Cxf for any function / with support in K where Ck is a 
compact operator in L2 (K) . 

Proof of Theorem. 

Lemma 3.2 The integral transform 

Inf (x) = {P) f [ ^^J^ ^1 dAdy?, for even n 



n — 1 J\ 

n-1 



Inf{x)= I ^^Mf{X,u] 



is a continuous operator L2 (-^)comp ~^ -^2 (-^)ioc • 

A proof of the Lemma is given in P^. For even n and an arbitrary x G X and a 
function / that vanishes in a neighborhood of x we calculate 

Inf X = P / — — — - — ^ — = / dw P / / / [y) q— ^ 

is <^"'(x;A,a;) Jgn-i JmJe(y,u^)=\ (^(x,a;)-A) 



(P) I 7777 ^1 fdeAq= [ Re e„ (x, y) f (y) dV (y) 
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Here, the relation dA = dO holds in Z and the equation dO A q = dV is fulfilled in X by 
definition. Thus the function 0„ is the off-diagonal kernel of the operator /„. It vanishes 
since of the assumption. Therefore B„ {x, y) is supported in the diagonal and according 



to Lemma [372] we have 6„ {x, y) = an (x) 6x (y) for a locally bounded function a„ in X. 

{e {x, u) - A) M$/ (A, u) dXdu 



If n is odd we have 
U (x) = 



du 



d 



in— 1 



(9A"- 



f (y) Q = I 

zix,u) J X - e{y,u) - tO 



fq\ = I Qn{x,y)fiy)dV{y) 

Z(X,uj) y X=0(x,ui) 

f {y) dV 



d 



in— 1 



dX 



n-l 



fq 



dV 



{n-l)\ 

7T 

(n-l)! 
7r 



Im 



9{x,co)-9{y,uj)-i0f 
Ime„ {x,y) fdV 



du fdV 



that is vr^^ (n — l)!ImO„ is the kernel of J„ which vanishes off the diagonal by the 
assumption. By Lemma 3.2 we again conclude that {n — 1)! Im9„ = a„ {x) 5^ iv) for 
a locally bounded function a„. 

Next we calculate the function a„. Choose a smooth function cq of one variable with 
support in [—1, 1] such that cq (0) = 1 and set (x) = cq (|a;|^ /e^) for x G MJ^ and any 
e > 0. Take a point € X and show that 



Re / dV 



Ce (x — Xq) dcj 



for n even and 



(n-l)! 



Im / dl^ 



(x, uj) — 6 (xo, w) — iO) 
(x — Xq) duj 



Co (|x| /e^) for x G 



-> ctn (a^o) 



^ Jx Js"-^ {(^ ix,Lu) — 6 {xo,uj) — iO) 

for n odd as e — )■ 0. We can change order of integrals and integrate first over X. 
Lemma 3.3 // n is even, we have for any xq G X, arbitrary u G S"~^ and small e 

e, (X - Xo) dV (_l)"/2-l T,(n+l)/2 ^ 



an {xq.uj) = Re 



X 



(x, u) — 9 (xo, co) — iOy 



r((n + l)/2) \V9{xo,cu) 



where o(l) < Ce^^'^logl/e where C does not depend on u. For odd n we have 



an [,Xq,UJ) 



[n 



IV. 



Im 



TT 



Ce (x - Xo) dV 



X 



'{x,u) -9{xo,co) -iO)" (n-l)! \V9{xo,uj] 



+ 0(1) 
(4) 



Tn+0 (e) 
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Taking the limit and integrating ^ over S" ^ yields for even n the equation 

an (xo) = lim [ a„ {xo,u) du = (-l)"/^-! J f 

e^oj T{{n + 1) /2) J |V6'(xo,a;)| 

(2.,)" 1 /■ d. _ 



(n- 1)! |S"-i| J |V^(xo,w)r (n- 1)! 
which implies (|2]). For odd we obtain 



a„ (xo) = lim / a„ (xq, dw = 2 {2ni)" ^ — ^ /" — ^ = 2 (27ri)" ^ L>„ (xq) 



This yields ([3]). This completes the proof of Theorem 3.1 ► 

Proof of Lemma. We show first that the 6 can be replaced by a linear function. 
Assume for simphcity that xq = and d9{xo,uj) /dxi = \Vg9\ ,d9{xo,u) /dx2 = ... = 
d9 (xo,a;) /dxn = for a coordinate system xi, ...,x„ in a neighborhood of Xq such that 
g (xo; d/dxi, d/dxj) = 6ij. We have then dV = vdx, dx = dxi...dx„ where f is a smooth 
function such that v (xq) = 1. 

If n is even, we integrate by parts n times with respect to Xi : 

f (x - Xq) dV ^ 1 r di (x) dx 

y (^(x,a;)-^(xo,c^)-^Or (n - 1) J {9 {x,co) - 9 {xo,uj) - lOf-^ ^' 



1 



log \9 {x,u) — 9 (xo, w) I dn (x) dx 



di 



d (x — Xo) f (x) dce V {di\Vg9 {x,u!)\ — vdi\Vg9 {x,co)\) 



dxi \Vg9{x)\ dxi\Vg9{x,io)\ |Vg^(x,a;)| 
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^ d dn-i (x) ^ ^"6^ 

" ~ 9xi |Vg^(x)| ~ 9x^ |Vg^(x,w)r 

where omitted terms only include derivatives of of order < n. Changing the variables 

X = ey we get 

dn (y) = + o (.^-") (6) 

dy'l\Vg9{ey)\ 

By Lagrange's theorem we can write 

9 {ey, u) - 9 (0, u;) = ep, (y) , (y) = I {y, V9 {sty)) dt (7) 

Jo 

This yields 

log \9{x,u) -9 (xo, io) I dn (x) = j (log ^) (a::) dx + j log jp^ (y) | (i„ (x) dx 
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The first integral vanishes since dn is equal to Xi-derivatives of a function with compact 
support. It follows that the left hand side of ([s]) equals 

j log \Ps (x) dx = J log Ip, (,)| |v.(,^,^)p g^dy + O is) 



{n-l)\J V— j 'iVg^^K,^)!"^^^ 

since the logarithmic factor is absolutely integrable. By ([7| we have C^-convergence 
— )■ (y, V6' (0)) as e — 7- in a neighborhood of the origin. This implies the inequality 



log \yi \dy- log \pe (y) \dy <C6 log \6\ 

l<i,|yi|<<5 Ay\<'i-My)<s 

where 6, < 5 < 1 is arbitrary and C does not depend on e and 6. On the other hand, 
log Ipe {y) \ — > log \yi\ everywhere as e — > 0. Therefore 



j log \pe {y)\ Q-n^y ^ j log Q-n^y 



and 



(^(x,a;)-^^(xo,u;)-^Or (n-l)!|Vg^(xo)r y l^^l Sy^^''^ 

More detailed arguments show that the difference is equal to O [e^/"^ log e) . The same is 
true for the linear function 6 (x, u) = X\ that is 

ee{x-xo)dV 1 f e (y) dy 



x'\Q{x,bj) -Q{xQ,bj) -i^Y |Vg6' (xo,c<;)r A" (z/i - ^0) 
Calculate the integral in the right hand side by partial integration: 



Re / , ''^y =^Re ^ 



{yi -iOf n-1 V (yi - iO)" ^ dyi 

— -Re f (coswi — iO)^~" dcj f ^~^dr^ 

Re / (cos Wi — zO)^~" do; 



n 
1 



n — 1 



n-l 



n — 1 



Re / sin"~^ cji (coswi — iO)^~" dwi (8) 
Jo 



where y = r cos ui since 



de 
dy 



- = 2yi—^, yi — iO = (coscji — iO) r, / deo/dr'^dr'^ = — e (0) = —1 
1 c/^ Jo 



By substituting s = cos^ Ui we get 



n-2 

I _ 

coswi + zOy 2 V 2 ' 2 / ' 2 



Re/ (^^^) " dc., = l5(^,^)=(-ir/-^" 



We use the formula for Beta-function extended for all complex (non negative integers) 
values of arguments. The exponent A = l/2 — n/2isa regular point and we can use a 
classical formula. The right hand side of ([8| equals 



2(n-l) ' ' r((n + l)/2)' 
which yields Q. 

In the case of odd n we integrate by parts as in the previous case and get 

Im/ jx - xo) dV = [ dn(x)dx 

J {0{x,uj) -e{xo,uj) -zO)'' - 1)! ye(x,<^)>e(xo,a;) 

Taking in account ^ and convergence — )■ yi the limit of the right hand side is 

1 / ^dy (9) 

Integrating by parts backward gives the equation 



where y' = (j/2, ■■■,yn) , m = (n — 1) /2, s = r^. Here, we applied the formula 

d^'-'eiy) =2-(n-2)!!-^'"'°^'^ 



yi=o 



dyr' 

Integrating by parts m — 1 times in the interior integral we get the quantity 



' '9'"eo (s) ^ 1 _ 1)! / 9eo/9sds = (-I)'" (m - 1)! 



This implies that ^ equals 

(_l)-2"^-i7r(m-l)!(n-2)!!|S"-2| i 2 (27ri)"-' 



(^-1)! \V,9{xo,oj)r {n-l)\ |S-i||Vg^(xo,a;)r 

For odd n we have 

^ (^-1)!k^t^ /- e,(x-xo)dV _ 2(2vr.r^ 

{e{x,oj) -e{xo,u) -tof |S"-i| |Vg^(xo,w)r 

2(27rz)""^ f du 



a„(xo,w) = ^limlm f 

Integrating over S"~^ we get 



and ([3]) follows. I 



4 Integrals of rational trigonometric functions 



We focus now on the condition Q {x,y) =0. A function of the form 

k 

t (ip) = aj cos + bj sinjip 

j=0 

is called trigonometric polynomial of degree k if ak ^ or bk ^ 0. Any trigonometric 
polynomial is 27r-periodic and is well-defined and holomorphic on the cylinder C/27rZ. 
It always has 2k zeros in the cylinder. If the polynomial is real, number of real zeros is 
even. 

Lemma 4.1 Let t{(f) and s {(p) be real trigonometric polynomials such that deg s < 
deg t and all the roots oft are real. Then for r = s/t and arbitrary natural n 



(P) r" (y,) = 2 (^) + + 2 X ~ ^ °' 

Proof. Suppose first that all roots of t are simple. Let ai < 0:2 < ... < oim be all 
roots of dt/dif on the circle ]R/27rZ. Let = sgn dt/d(p on the interval (0^,0^+1) for 
k = 1, ...,m, where ttm+i = «!• The function r (C) is meromorphic for ( = ip+ir G C/27rZ 
and has no poles in the half-cylinder {r > 0} since of the assumption. Take a continuous 
function A = r {(p>) defined on the circle that vanishes in all points and is positive 
otherwise. We have for any k 

Qk+i rotk+1 

{r {(p) + EkiOT dip = r" ((^ + iO) d(^. 

Ofc ''oik 

Take the sum 



I (r (^) + Ski^f dip = J2 [ ^" (V' + ^0) = / ^" (V' + ^0) d^ (11) 

^Ofc Jo 



Now we can replace the form r" + iO) dy? by r" ((^) dC for C, = ip + irj ioi an arbitrary 
> without changing the integral in the right-hand side. We have — )■ as 
— > 00 hence the right hand side of (11) vanishes. The real part of the left hand side is 
equal to the left hand side of (10) which also vanishes. 

In the general case, we can approximate the polynomial t with real roots by polyno- 
mials t with real simple roots. The equation (10) holds for f = s/t hence it is true for 
r = s/t. ► 

Proposition 4.2 Let v G M" and a eM. be such that \a\ < \v \ . Then for arbitrary even 
n>2 

Re / — ^ = 

{{00, V) — a — tO) 
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and for arbitrary odd n > 3 



Im 



0. 



{{co,v) — a — iOy 
Proof. We may assume that |f | = 1. For even n we have 



Re 



do; 



{{co,v) — a — iO)^ 



Re 



do; 



(cos (f — a — iOy 



where (p is the spherical distance between cu and v. We have do; = sm'^~'^ ipdipdu' where 
du' is the area of a unit sphere S"~^ = {w e S"~^; ip = vr/2} . Integrating over n — 2- 
spheres (p = const we get 



Re 



du 



2n-2\ 



(cos (f — a — iOy 



2tt „• n-2 



sin" {pdip 
(cosy^ — a)" 



since the integrand is vr-periodic. The right hand side vanishes by Lemma 4.1 
For odd n we have 



Im 



du 



2n-2| 



(cosv9 — a — iO)" 2i {n — ly. 



sin" Lfdip 



sin" Lfdif 



Jq (cos v9 — a — iO)" (cos (y9 — a + iO)" 



sin" ^ ipd(p IT |S 



n-2| 



sin" ^ ipdip 



res. 



2z(ra- 1)! (cos(^-a)" (ra - 1)! (cos(^-a)" 

where a = arccosa G [0, tt] . Changing variable ( = cos if and omitting the constant 
coefficient we get the quantity 

(1-CTdC 
"'"''^ iC-af 

where m = (n — 3) /2. The residue is equal to zero since the numerator has order 2m = 
n-3. ► 



Corollary 4.3 For any regular resolved generating function ^ = 6 — X such that for any 
pair of points x y in an open set Q G X we have 9 (x, uj)—6 {y, u) = {v, u)+a, \a\ < \v\ , 
formulas M) and hold for any function f E L2 (X) with support in ^l. 



5 Reconstruction in spaces of constant curvature 

We apply the above results to recover few known and unknown inversion formulas for 
geodesic integral transforms in spaces of curvature k = 0, 1, — 1. 

Euclidean space. Take the generating function $ (x; A, u) = {u, x) — A in M" x 
E, E = M X S"-^ We have |V^| = Dn (x) = 1. Then (g,^ coincide with the classical 
reconstruction in Euclidean space from data of hyperplane integrals. 
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Elliptic space. Funk |[5] inspired by tlie seminal paper of Minkowski [12j found a 
reconstruction formula of an even function / on the unit sphere from integrals Gf 
for the family of big circles. His method was adapted by Radon for Euclidean plane. A 
generalization of Funk's formula for arbitrary dimension is due to Helgason [8]. 



A reconstruction formula of a different form can be obtained if we apply Theorem 3.1 
to a generating function $ (x; A, u) = {u, ~ ^ defined in X x S where X = {{xq, x) G 
E"+\xo + \xf = l,Xo > 0}, S = G -E" : 1^1 < 1} and g is the metric in X induced 
form the Euclidean space E"''^^. Omitting some simple calculations, we come to 

Theorem 5.1 If n is even then any function / G L2 (-'^)comp reconstructed from 

its integrals Gf (a) over big spheres S (a) = {x E X , {a, x) = 0} , cr G S" by 

(n-1)! f Gf{a)da 



(2vri)" Js^ {(T,xy 
where S"~^ = {a G IR"'^^; \a\ = l,o"o ^ 0} is a hemisphere. If n is odd we have 



f (x) = — / ^^""'^ iW^ ^)) Gf (a) da. 



2 {2mf-^ 



n-1 



Hyperbolic space (k = — 1) • Take the generating function $ (x; A, w) = — A, 6 = 
—2 (jxf + 1) (w, x) , — 1 < A < 1 in the unit ball X C M". The hypersurfaces Z {X,uj) 

(1 2 \ — 
1 — |x| ) ds. By a similar rear- 
rangement we obtain 

,(.)._(!Lzi^/ » (12) 

i^^V Jq+ 

for even n and 

/ (x) = ^ / 5^"-'^ {{cr, x)) Gf {a) da (13) 

2 {2Tn) Jq+ 

for odd n where = {cr G ]R"+^; a^ — al — ... — a'^ = —\,aQ>Q} is the dual one-fold 
hyperboloid. 

A reconstruction in a different form was done first by Radon p!7] {n = 2) and Helgason 



[8] (n > 2); formulas (12) and (13) are due to Gelfand-Graev-Vilenkin |7j- 



6 Equidistant spheres and horospheres in hyperbolic 
space 

Equidistant spheres. Let X be again a unit n dimensional ball, n >2 and 

$ {x; X,uj) = e-X, e {x, 00) = ^~ ^^'f , CO G s"-^ 

1 — 
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be a generating function where < p < 1. For a fixed oj and an arbitrary A 7^ the 
hypersurface Z (A, u) = {x; $ (x; A, u) = 0} is the intersection of X and of an n—1 sphere 
S (A) whereas S (0) = {{u , x) — p = 0} is a hyperplane; all the spheres S (A) contain n—2 
sphere S (0) n dX. For arbitrary real A and /i the hypersurfaces S (A) nX, S (n) r\X are 
equidistant with respect to the hyperbolic metric. Check that $ fulfils the conditions of 



Corollary 4.3 for arbitrary p, < p < 1. A proof of regularity is a routine. Further we 
have 



e{x,uj) -e{y,uj) = - {uj, 



and we need to prove that 



X 



\x\ 



1 



x\ 



X 



1 - \y\ 



> 



(14) 



for arbitrary x ^ y E X. Squaring both sides we reduce (14) to the obvious inequalit'^ 
2{l — xy) > (2 — |a;|^ — The proof is complete and reconstructions (2) and (3 

follow. Write these formulas in terms of the geodesic integral transform 



Hf{a)= [ fd,S, a={X,u). 

Jz(a) 



where dgS is the hyperbolic hypersurface element. The Minkowski- Funk operator M$ 
can be written in terms of the Euclidean integral transform Gf (cr) = /^(-^-j fdeS since 

of the factorization |V^(x,w)| = (1- |a;|^) v^4A2 - ApX + 1 (see ill). On the other 
hand, 

2 



which yields 



where 



Mf{\,u) 



,1 — 
Gfi{X,u) 



n-l 



d.S 



v/A2-pA + l/4 v^A2 -p\ + 1/4 



il - \x?^ 
fi (^) = ^ J V (x) , h {x) 



1 — \x\ 



fix). 



Corollary 6.1 For any function f with compact support in the unit ball, a reconstruction 
is given for even n by 



fix) 



1 {l-\xf) 



2\2n 



XI 



Hf{\,uj) d\ 



du 



v/A2-pA + l/4 {{co,x) - A (1 



1x1^))" 



(15) 
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and for odd n by 



1 (l-l^l')'" 



d 



n-l 



Hf{\,L0) 



2 (47ri)" (x) Jsn-i 9A"-i -pA + 1/4 



do;. 



(16) 



X=e{x,u;) 

Horospheres. Taking p = 1 in the above formulas we obtain the function 

1 — {u, x) 



9 (x, uj) 



\x\ 



which defines the family of horospheres 6{x,uj) = A, 1/2 < A < oo. The formulas (15) 



and (16) holds for horospheres if we substitute y/AX^ — 4A + 1 = 2A — 1 and integrate in 



(15) over the ray (1/2, oo) in the interior integral. 

For the cases n = 2 and n = 3 reconstruction formulas (of different form) are contained 
in the book of Gelfand-Gindikin-Graev [6j. 



7 Isofocal hyperboloids 



Hyperboloids. The equation A = |x| e > 1 defines a fold of a two-fold hyperboloid 

2 ' 



^AA 2 2 

axi — Xo — ■■■ — X = a 

a J a"' 



with a focus at the origin. The function $ (x; A, w) = 6* (x, u) — \, 9 (x, ip) = \x\+e (x, u) 
generates the family of all one-fold hyperboloids in X = M" \ with a focus at the origin. 
The function 

9 (x, u) — 9 {y, u) = e {u,x — y) + \x\ — \y\ , u E S"""^ 



satisfies the conditions of Proposition 4.2 since ||x| — \y\\ < e\x — y\ for any x,y E 



X ^ y. Therefore Theorem 3A holds for this family. We have 

\V9 {x,uj)f = l + e^ + 2e\x\~^ {u,x) 

and 



S-n-i (1 + ^2 + 25 



X {U,X 



where (p is the angle between u and |x| ^ x. This integral does not depend on x. 

Corollary 7.1 For any smooth function f with compact support in M" and even n, the 
equation holds 

1 r f Mf{X,oo)dXdoo 



)) 



n/2 



{1 + e'^ — 26 cos v?) 



n/2 



fix) 
fix) 



(27rz)"L'„ 7g 
1 

2 (27ri)"~^ Dn 7s"-i "^A 



(|x| + e {uj,x) - A) 

Qn-l 



-,Mf{X,uj)\,__ 



n even 



|x|+e(x,<^> dw, n odd. 
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8 Photoacoustic geometries 



Consider a resolved generating function 

$ (x; A, = |x - e {uj)f - A, w G S"-i 

where ^ : S""^ — )■ M" is a smooth map. We call the image C of ^ central set. Any surface 
Z (A, 00) is a sphere of radius with the center ^ (uj) G C and 

M,/(A,c.)= 



2yA 

where G/ (A, u) is the Euclidean integral over this sphere. Inversion of the operator 
M$ implies inversion of the spherical integral transform G for the given central surface 
C (and vice versa). This subject is of special interest in view of application to the 
photoacoustic (thermoacoustic) tomography, see surveys [IH],[n]- Inversion formulas for 
a function supported in a half space with the hyperplane central set is found by Fawcett 
P]. In [ll] a reconstruction was done by reduction to the Radon transform. For a 
spherical central surface Finch with coauthors [3],|1] found a reconstruction formula of 



type (17)-(18) in the physical domain for arbitrary dimension. Another reconstruction 



formula was proposed by Kunyanski [lOj; it is similar to (17)-(18) after a simplification. 
An inversion of the spherical mean operator in three dimensional space was constructed 
by Xu and Wang [19] in terms of a Dirichlet-Green function. This approach gives an 
explicit reconstruction in the frequency domain when C is a sphere or a circular cylinder. 
Kunyanski [11] construct inversion for polyhedral center sets with special symmetries. 
Natterer [TB] has found an explicit inversion in the physical domain when the central 
surface C is an arbitrary ellipsoid. Surveys of related results are done in [1] , [19J , [Hj . We 
show here that for an arbitrary ellipsoid (elliptical cylinder) in i?" a simple reconstruction 



is given by Theorem 3.1 



Ellipsoids. Set Ci<^) = {o-i^i, ■■■,ctnU!n) where positive constants. The 

central hypersurface {x = {00) ,uj G S"~^} is the boundary of an ellipsoid Ea with half 
axes ai, a„. Then 

e {x, u)-e {y, u) =2{^{u) ,y - x) + - \yf = 2 {u, z) + \x\^ - \yf 

where z = (ai {yi — Xi) , ...,an {ijn — a;„)) . The inequality holds 

- \y?\ = E " iyi + ^i)\ < X] l«i ivi - h*"^ {yi + Xi)\ < \\z\\ \\w\\ 



where w = [a^^ {xi + yi) , [xn + ?/„)) . Suppose that x,y G Ea and x ^ y then the 
point [x + y) /2 belongs to the interior of Ea which implies ||w|| < 2. It follows that the 



right hand side of is strictly bounded by 2 \\z\\. By Proposition 4.2 Theorem 3.1 holds 
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for any n > 2. Note that |V6'| = 2 |a; — ^ (a;)| = 2y/X. It follows that any function / 
supported in the closed ellipsoid Ea can be reconstructed by the formula 

f (a;) = ^ f (P) f ^^(^''^^^^^^ (17) 

for even n where we did the substitution A = p^, and by 

1 f fdy-'Gfip^u) 



for odd n where 



do; (18) 



^"^"'^ 2"is-i|/ ix-eMr ^^^^ 

Cylinders. If a central set C is unbounded can not apply the same method since it 
can not be regularly parametrized. However one can write a reconstruction formula for 
any closed cylinder E with an elliptic base. Indeed, E is a union of the family of ellipsoids 
Ea as one or several half axes, say ai, Op, tend to infinity, ap+i, a„ being fixed. One 
can come to limits in (17) and in (18). We omit details. 

Algebraic plane curves. In the case n = 2, there are more geometries which admit 
exact reconstruction formulas. We call a curve C C trigonometric of degree k if it is 
given by a parametric equation 

where ^i,.^2 are real trigonometric polynomials of degree k. A trigonometric curve is 
always a component of a real algebraic curve. A point x G is called hyperbolic with 
respect to a trigonometric curve C if any straight line L through x meets the curve at 2k 
points (counting with multiplicities). It is easy to see that the set H (called hyperbolic 
set) of all hyperbolic points is always closed and convex. Introduce an Euclidean structure 
in and consider a function 

9{x,^) = \x-av)\\ e(^) = (6(<^),6(<^)), ¥^eS\ xeR^. 

Proposition 8.1 Let H be the set of hyperbolic points with respect to a trigonometric 
curve of degree k. For arbitrary points x,y E H, x ^ y, all roots of the polynomial 
6 (x, (p) — 9 {y, ip) ( of order k ) are real. 

Proof. We have 

e{x,^)-e{y,(^) = \x - ^ {ip)f - \y - ^ {ip)f = \xf - \yf -2{x-y,^{ip)) 

= 2{x-y,av)-s) (20) 
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where s = {x + y) /2. This point is contained in H since H is convex. Therefore, the hne 
L = {z = s + rv,r G M} has 2k common points (v^i) , (v^2fc) with C for arbitrary 
vector f 7^ 0. If f is orthogonal to x — y the right hand side of ( [20| vanishes. The 
corresponding angles v^i, '^2k are real roots of the polynomial 9 {x,ip) — 9 {y,ijj). ► 
The family of circles centered at the curve C is generated by the function $ (x; A, ip) = 



9 {x, (p) — X. Applying Lemma 4.1 we get 



Corollary 8.2 Let C be a trigonometric curve with a hyperbolic set H. Theorem 3.1 
holds for the family of circles centered at C, arbitrary function f supported in H and for 
both integral transforms M$ and G. 

There is a large variety of trigonometric curves C with non empty hyperbolic sets. 
Examples. 1. Let [ip) = 2 cos 2(f — cos (f, ^2 (v^) = 2 sin 2(p + sin (p. The curve C is 
shown in Fig.l 



Fig.l 

The hyperbolic set H is the triangle in the middle. 

2. A hyperbohc "square" set is defined by the trigonometric curve ^1 (ip) = 2cos3(y9 + 
COSV9, ^2 if) = 2sin3(y9 — simp, see Fig.2 
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Fig.2 



3. A "pentagon" is the hyperbolic set of the curve {(p) = 5 cos 4(^+4 cos ip, ^2 {f) 
5 sin 4</7 — 4 sin (/?, see Fig. 3: 
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